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Abstract. We study the impact of different cold nuclear matter effects both on J/ip and ip' 
production, among them the modification of the gluon distribution in bound nucleons, commonly 
known as gluon shadowing, and the survival probability for a bound state to escape the nucleus 
-the nuclear absorption. Less conventional effects such as saturation and fractional energy loss 
are also discussed. We pay a particular attention to the recent PHENIX preliminary data on 
ip' production in dAu collisions at y/s = 200 GeV, which show a strong suppression for central 
collisions, 5 times larger than the one obtained for J/ip production at the same energy. We 
conclude that none of the abovementioned mechanisms can explain this experimental result. 



1. Introduction 

The properties of the production and the absorption of quarkonium in proton(deuteron)-nucleus 
collisions, such as J/ijj and suppression in dAu collisions at RHIC, offer the opportunity to 
study Cold Nuclear Matter (CNM) effects in the nuclear medium. This is of particular relevance 
for the understanding of QCD at high density and temperature in nucleus-nucleus collisions, since 
it is in fact adventurous to claim that the charmonium family is a quark-gluon plasma (QGP) 
thermometer [1] without calibrating first the CNM suppression. 

We have studied the impact of different cold nuclear matter effects both on J/r/j and ip f 
production, namely the survival probability for a bound state to escape the nucleus -the nuclear 
absorption- and the modification of the gluon distribution in bound nucleons -gluon shadowing. 
Alternative mechanisms such as saturation and fractional energy loss are also considered. 

We have then discussed whether these CNM effects could explain (or not) the data on J/ip 
and ip' within a coherent picture. We have paid a special attention to the recent PHENIX 
preliminary data [2J on production in dAu collisions at \/s = 200 GeV, which show a strong 
suppression for central collisions in the mid-rapidity region. 



2. Propagation in the cold nuclear matter: J/ij) vs tf 

The probability for the heavy-quark pair to be broken up in the propagation through the nuclear 
medium is commonly known as the nuclear absorption. It is usually parametrised by an effective 
break -up cross section cr e ff. 

In principle, the larger size of the tf compared to the J/ij) tells us that the tf should suffer 
more break-up than the J/ij). Nevertheless, it is important to remember that the relevant 
timescale to analyse the pair evolution is its formation time, which, following the uncertainty 
principle, is related to the time needed - in their rest frame - to distinguish the energy levels of 
the IS and 2S states, t f = (m2 2Mc ^ 2 ) = 2 x 3.3 GeV / 4 GeV 2 = 0.35 fm for the tp. Moreover, 

tf has to be considered in the rest frame of the target nucleus, i.e. the Au beam at RHIC. 
The boost factor 7 is obtained from the rapidity of the pair corrected by the Au beam rapidity 
7 = cosh(y - y£u j where y Au^ = _ 5M for RHIC> 

The formation time for the different rapidities reached at RHIC are given in Table [[} As it 
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Table 1. Boost and formation-time ^-dependence in the Au rest frame of the ip at ^/snn = 200 GeV. 

is shown there, in the mid and forward rapidity regions at RHIC, tf is significantly larger than 
the Au radius. This means that the cc is nearly always in a pre-resonant state when traversing 
the nuclear matter. This implies that the break-up probability should be the same for the tf 
and J/^, since these states cannot be distinguished at the time they traverse the nucleus. 

3. Gluon shadowing and antishadowing 

At high energy, the nucleons in a nucleus shadow each other, leading to a modification of the 
nuclear Parton Distribution Functions (nPDF). This modification, in the form of a shadowing- 
correction factor, has to be applied to the J/ifj or tf yield in dAu collisions. This shadowing 
factor can be expressed in terms of the ratios Rf of the nPDF in a nucleon of a nucleus A to 
the PDF in the free nucleon. 

We shall consider three different shadowing parametrisations: EKS98 [3J, EPS08 [4J and 
nDSg [5] at LO, which have been implemented in our Glauber Monte-Carlo framework JIN [6J, 
already used to describe J ftp production at RHIC [3 El [9]. Our code accounts for improved 
kinematics corresponding to a 2 —> 2 (g + g — >> cc + g) partonic process for the \j) production as 
in the Colour-Singlet Model at LO [10] or the Colour-Octet Mechanism at NLO. 

In Fig. [TJ the effect of the nPDFs on J/ip is compared to PHENIX data on J/ip and tf 
production [2J. We have used the three shadowing parametrisations mentioned above, together 
with four possibilities for the effective nuclear absortion, namely <r e ff — 0, 2, 4, 6 mb. 

Note that, while the choice of identical break-up cross sections for J/i/j and tf is perfectly 
justified given the argument of large formation times invoked above, the common prejudice 
that there is a similar shadowing in both cases needs further explanations. It is true that the 
shadowing corrections, that depends on x and /ip, should be very similar for both particle 
production, due to their similar masses. Nevertheless, the difference between their masses -of 
the order of 0.6 GeV- could induce both a shift on the corresponding x at a given rapidity and 
on the corresponding /ip. We have taken the factorisation scale, ^p, equal to the particle mass, 
3.7 GeV for the tf and 3.1 GeV for the J/ip. In Fig. § our results for the ratio between the 

shadowing factors — jf^ for different shadowing parametrisations -EKS98, EPS08, nDSg at LO 




Figure 1. Nuclear modification factor in dAu collisions, RdAu, for J/ip and ip' as measured by PHENIX 
at <s/snn = 200 GeV versus the number of collisions for the mid-rapidity range \y\ < 0.35 compared to 
our computations for four values of the nuclear absorption (from top to bottom: a e ^ = 0, 2, 4, 6 mb) 
using EKS98 (left), EPS08 (middle) and nDSg (right). 



along with EPS09 at LO and NLO pT] and nDSg at NLO- are shown. Note that the difference 




Figure 2. ip' over J ftp nuclear modification factors in the rapidity range —2.2 < y < 2.2 (left) and in 
the mid-rapidity region \y\ < 0.35 (right) using EKS98 (black), EPS08 (red), nDSg (green) and EPS09 
(blue) at LO (continous line) and NLO (discontinous line). The shaded zone indicates that tf is there of 
the order of the Au radius. In this region, the J/ip and ip' absorption may be slightly different. 



of the shadowing impact is at most 5 %. Even more important, this difference will favour ip f 
production in the mid-rapidity region, in contradiction to what experimental data show. 

4. Other effects: gluon saturation and fractional energy loss 

The possibility of a gluon saturation as predicted in the framework of the Color Glass Condensate 
(CGC) \12\ [T3] is sometimes invoked to have an impact on quarkonium production. We have 
already discussed [H] that this impact would be negligible in the case of T production at RHIC 
energies. Let us now discuss if it could affect J/ip and ip f production. The relevant scale below 
which one expects significant non-linear effects in the evolution of the gluon distributions in pA is 

the saturation momentum. It can be evaluated as [15] Q 2 sA = is x0.2x (in units of GeV 2 ), 

with A rsj 0.2 -T- 0.3 and xq = 0.01. Typical values for RHIC kinematics are given in Table. [2] 
together with the ratio of Q s to the ip mass. Note that in the backward rapidity region x is always 
above xo -the maximum momentum fraction below which parton saturation can be considered. 
In the mid and forward rapidity regions, the saturation scale is always below the typical energy 
scale of the process, namely m^f or mj/^. In particular, in the mid-rapidity region the ratio of 
Q s to the ip mass is |, which makes unlikely the possibility of specific saturation effects. Only 
in the most forward rapidity region there is a small probability for phenomena beyond collinear 
factorisation i.e. beyond those encoded in the nPDFs. Moreover, these effects -if they exist- 
would suppress the J/tp slightly more than the ip\ in contradiction to the experimental results. 

Recently, the possibility of a medium-induced energy loss that scales with the energy E has 
been pointed out [16] . This fractional energy loss would affect the heavy-quark pair produced 
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Table 2. Evaluation of the saturation scale in the Au nucleus, Q S Au> a t \/snn = 200 GeV in dAu 
collisions as a function of the ijj' and J/ip rapidities. 

in a coloured state provided that it has long enough time to radiate. The fraction of medium- 
induced radiated energy is given by AE/E = Ax\/x\ ~ N c a s A (p^ ) / Mt , where A(pj) is the 
broadening of the radiated gluon from the proton and Mt is the transverse mass of the final-state 
coloured object. If 100% of the ^'s were produced in a coloured state at RHIC, Ax\/x\ ~ 10%. 
Assuming that the corresponding suppression factor is just the ratio between the PDF of the 
gluon evaluated at its original x' 1: x[ = x\ + Ax\ before the loss and the PDF at the resulting 
Rloss(%i, Q 2 ) — g( x i, Q 2 )/ d{ x ii Q 2 )- Taking modern global- fit gluon PDFs, we find that this 
fractional energy loss could give an extra suppression of the order of 10% in the mid-rapidity 
region and up to 15% in the forward region. Note that this suppression only applies when the 
-0 is produced in a coloured state, and that its effect is similar in J/0 and 7 production due to 
their similar masse^J Thus, this cannot explain the experimental found either. 

In conclusion, we have studied different CNM effects and their impact on J/0 and 7 
production. We have shown that, at RHIC energies and in the mid-rapidity region, the results 
from the four CNM effects which we investigated, namely an effective nuclear absorption, the 
gluon shadowing, the gluon saturation and a fractional energy loss, impact similarly the J/ip 
and the 0/, in contradiction to the recent PHENIX preliminary data on 0/ production in 
dAu collisions at 200 GeV. The abovementioned effects should be taken into account in the 
interpretation of the upcoming LHC data on J/0 and 0/ production. 
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